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We present a theory of the influence of band renormahzation and excitonic electron-electron in- 
teraction effects on the optical conductivity oiua) of doped bilayer graphene. Using the Keldysh 
formalism, we derive a kinetic equation from which we extract numerical and approximate analytic 
results for o"(w). Our calculations reveal a previously unrecognized mechanism which couples the 
Drude and interband response and renormalizes the plasmon frequency, and suggest that screen- 
ing must play an essential role in explaining the weakly renormalized conductivity seen in recent 
experiments. 



PACS numbers: 73.63.Bd,78.67.Pt,72.20.Dp 

Introduction — Experimental progress flj over the past 
five years has made it possible to isolate single-layer 
graphene (SLG), an atomically two-dimensional electron 
system based on a honeycomb lattice of carbon atoms, 
and study its electronic properties. More recently similar 
techniques have been used to study multi-layer graphene 
films and other types of graphitic nanostructures. One 
of the surprises in this field is that the electronic proper- 
ties of Bernal stacked bilayer graphene (BLG) are quite 
distinct [2| from those of SLG. The low-energy electronic 
excitations of a bilayer are massive and have momentum- 
space Berry phase 2tt, while those of SLG are massless 
and have Berry phase tt. 

The optical conductivity <t{uj) of BLG has received a 
lot of attention, both experimentally and theoretically. 
Theoretical studies ^ ^ ^ |3, 0] of (j{uj) in BLG have 
so far entirely neglected the interactions effects that are 
known to be crucial in the optical response of regular 
semiconductors [1, 0|- The non-interacting electron in- 
terpretation of <j{uj) data [l^, 11, 12, 13| does never- 
theless appear to be generally successful, surprisingly so 
since BLG is generally expected to display stronger in- 
teraction effects than SLG because of its parabolic band 
dispersion. For example band structure renormalizations 
are expected [l3| to be relatively modest in SLG, but 
substantially stronger in the BLG case. Indeed the in- 
teracting electron problem in bilayer graphene poses a 
number of interesting new questions because of its unique 
massive chiral quasiparticles. 

These circumstances call for the theoretical analysis 
of the influence of interactions on a{uj) presented in this 
Letter. We use a quantum kinetic equation (QKE) de- 
rived using the Keldysh formalism to take e-e band renor- 
mahzation and excitonic effects into account on an equal 
footing (thus correctly guaranteeing gauge invariance). 
We show that although the energy dispersion of BLG is 
parabolic, its optical properties are very different from 
those of regular semiconductors or semiconductor bilay- 
ers: i) A new coupling between the Drude [i.e., in- 
traband) and interband optical transition channels fol- 



lows from the chirality of the BLG band eigenstates; ii) 
the Drude-interband coupling (DIG) is responsible for 
a renormahzation of the leading order long-wavelength 
plasmon dispersion; and Hi) because of the the chirality 
structure, screening is responsible for an especially strong 
suppression of interaction effects. 

Quantum Kinetic Equation — States near the Fermi level 
of bilayer graphene are described by the two-band en- 
velope function Hamiltonian \^ H — — e^cr • n, where 
Efc = fc^/2m, n — (cos20fc, sin2(/>fc), and cr is the Pauli 
matrix vector which acts on layer pseudospin degrees of 
freedom. [We set h = 1 throughout restoring it only 
in the final expressions for (^{lo)]. This Hamiltonian is 



valid when vpk <C 71 where vp z 
particle velocity of SLG and 71 
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interlayer hopping amplitude. (We neglect trigonal warp- 
ing, which is important only at low densities and ener- 
gies.) Although conduction and valence band eigenener- 
gies have the same quadratic dispersion in regular semi- 
conductors and BLG, the eigenfunction properties differ. 
In the BLG case the conduction and valence band eigen- 
states are both linear combinations of 7r-orbitals, whereas 
in the regular semiconductor case the two bands have or- 
bitals with different atomic character. We will see that 
this property alone profoundly alters the cr(w) theory. 
Furthermore bilayer graphene is gapless (ijj . 

To incorporate band renormahzation and excitonic ef- 
fects on an equal footing, we derive a quantum kinetic 
equation for bilayer graphene using the Keldysh formal- 
ism and a first order exchange- interaction approximation 
for the interaction self-energy. Importantly the inter- 
action term in the envelope function is diagonal 



pseudospin labels at each interaction vertex. To obtain a 
kinetic equation, it is customary to employ a Wigner rep- 
resentation in which the relative coordinates r = ri — 
and T = ti ~ t2 in the Keldysh Green function [l8| 
are Fourier-transformed to obtain momentum and en- 
ergy variables k and e, and then perform a gradient ex- 
pansions with respect to the 'center-of-mass' coordinates 
R = {n + r2)/2 and t = (ti + 12) 12. The 2x2 dis- 
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tribution function fk is obtained by integrating the 
Keldysh Green function over energy. For the case of the 
bilayer graphene Hamiltonian we find that 



dfk 
dt 



e5-^+z[-efcn-o- + S]fc,/fc] =0, (1) 



where 



k' 



^k-k'fk' 



(2) 



is the quasiparticle exchange self-energy. The property 
that the 2x2 self-energy matrix at one wavevector is sim- 
ply an interaction weighted average of distribution func- 
tion matrices at different wavevectors is a consequence of 
the model's pseudospin independent interactions. 



We consider linear response to an ac electric field, 
S = 5oe-*"* and write fk = fi°^ + fi^^ where /f ^ = 
(1/2) ni?(l — ficr ■ n) is the equilibrium distribu- 

tion function. Here np is the Fermi function at zero 
temperature. Using Eq.([T|) we find that 



i efccr • n + S^,/, 



(1) 



Sk + ^y^ Vfc-fc' 



Jk' 'J k 



(3) 



where now refers to the self-energy evaluated using 
fk = /f ^ in Eq. @, and 



Sk = -{eS ■ fe/2) [^nF{ik^^/^k] (1 - /icr • n) + (1/fc) ^ ^inF{^k^.) (e£ x fe) • (cr x n) , (4) 

is the driving term of the QKE. The first term in Eq. Q drives intraband transitions and the second term interband 
transitions. The second term on the right-hand side of Eq. ([3|) accounts for changes in the self-energy in the non- 
equilibrium state. Because of this term, Eq. ([3]) is an integral equation which can only be solved numerically. The 
distribution function can also be expressed as a sum of intraband and interband contributions: 

(5-fc) [A{k) + a- ■ nB{k) 4- i (cr x n)^ G{k) + a,H{k)] + {Sxk), [i [a x n)^ C{k) + (j,D{k) + cr • nE{k) + F{k)] , 

(5) 

where the l,cr • n, (cr x n)^, and ct^ components of each contribution capture respectively changes in total density, 
conduction vs. valence band density difference, interlayer co herence, and layer polarization. 



Substituting Eq. ^ into Eq. ([3]) yields a set of eight 
coupled equations. We find that E, F, G and H are all 
identically zero, that A{k) = —B{k) = {ie/2uj)6{k — kp), 
and that C{k) and D{k) satisfy the following set of cou- 
pled integral equations: 

e 



ujC{k) + 5kD{k) = 6{k~kp) 



k 



5Y.'{k) 



(6) 



5kC{k)+ujD{k) = 6i(fc-fcF) [(5S'^'^(A:) + (5S'^'^(fc)](7) 

where 6k = 2e^. -f Tik+ — ^k- is the energy needed to 
create a vertical interband excitation, 

Y,Y.^k-k'0{kp - Xk'){l + ^iXcos2(Pk'k)/2, 



-'kfi 



X=± k' 



(8) 

is the equilibrium self-energy in band /i, (jjk'k — 4>k' — 4>k, 
and the non-equilibrium self-energy changes are 

JE^(fc) = ^Ffe_fe.cos0fc,fe7^(fc'), (9) 

k' 

5Y.^^\k) = ^Vl._fc'Cos0fc-fcCos20feaC(fc'), (10) 

k' 

,5S*''(fc) = -i^Ffc_fc-sin(/)fcasin20fc,feB(fc')-(ll) 



Eqs. (I6|)-([7|) are the equations of motion for the inter- 
layer coherence C(fc) and layer polarization D{k) com- 
ponents of the distribution function and describe preces- 
sion of valence-band pseudospins in effective magnetic 
fields due to the band-energy separation {5k) and to non- 
equilibrium self-energy corrections which favor interband 
coherence {5Ti'^'^) and layer polarization ((5S^). The sec- 
ond contribution to SYi'^ Eq. pTjl couples Drude and in- 
terband response (DIG). This DIG mechanism is one of 
the principle results of this Letter, ft appears because 
the Drude conduction band Fermi surface oscillation in 
an ac electric field changes the exchange potential expe- 
rienced by precessing valence-band pseudospins outside 
the Fermi surface. 



The current can be evaluated from the 
perturbed distribution function using J — 
5^ffseTr[X;fe(l/2){jfe,4^^}], where g^gs = 4 is the 
product of the valley and spin degeneracies and 
jf. = dH/dk = —{k/m)[{(7 ■ k) x — {cr x fc)^ y] is the 
current operator. It then follows from Eq. ^ that the 
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FIG. 1: (A) Renormalized Drude weight T) versus density 
n. The enhancement of the Drude weight increases with de- 
creasing density. (B) Band renormalization {T.k+ — '^k-)/£F 
a,t k = kp versus density. These results were evaluated with 
Coulombic electron-electron interactions and dielectric con- 
stant K = 1, corresponding to a suspended graphene sample. 



conductivity 

POO 

ct(w) = -(2e/7rm) / dkk'^ [B{k) + iC{k)] . (12) 

Renormalization of Drude Weight and Plasmon 
Frequency — Before discussing our general results we 
comment on the influence of DIG on the Drude weight. 
At low frequencies the non-interacting response is the 
out-of-phase oscillations of the Fermi surface with respect 
to the electric field which is captured by B{k) (x i/uj. 
Because of the self-energy correction JS"^'^, there is also 
an interband [C(fc)] response with the same frequency 
dependence [See Eq. ([7])]. When a momentum relaxation 
time T is added to the theory the i/to contribution to 
the conductivity evolves into a Drude peak contribution 
proportional to t/(1 — icor). The coefficient of this 
contribution is known as the Drude weight T). When 
DIG is included we find that for bilayer graphene 



FIG. 2: (Color online) Real part of the optical conductiv- 



{2e'^£F/T^h)f> T 
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(13) 



where the interaction induced Drude weight renormaliza- 
tion is given to leading order in by 



V = 1 



2'KmeF 



Here TZ{x) = {4:/l5x^){{x + l){x'^ - x^ + l)E[4a;/(x + 
1)2] - {x^ + l){x - l)^{x + l)K[Ax/{x + 1)2]} and IK,E 
are respectively complete elliptic integrals of the first and 
second kind. 

One important consequence of DIG is renormaliza- 
tion of the plasmon frequency. In regular semiconduc- 
tors with parabolic dispersion, the plasmon frequency 



ity Kea{ijj) 



frequency uj at density n = 10 cm for 



unscreened (black) and screened (blue/grey) e-e interaction. 
Results for screening with n = 1 (corresponding to suspended 
bilayer graphene) and k = 4 (corresponding to a bilayer on a 
Si02 substrate) are quantitatively similar and the blue lines 
show the results for n — 1. The disorder broadening is taken 
as F = 1/r = 0.1 meV (solid line) and ImeV (dot-dashed 
line). The ideal non-interacting case is plotted as the dashed 
thin red line. 



ujr, has no long wavelength interaction renormalization 
[l9| because of Galilean invariance. Since graphene sys- 
tems are not Galilean invariant, their plasmon frequen- 
cies are [2^ renormalized. Using the well-known relation 
between the optical conductivity and the polarizability, 
a{uj) — \iinq^Q[ie'^ujIl{q, to) / q'^], the real part of the po- 
larizability for vpq <^ Lu and lu ^ sp Ren(g, w) = 
{2£p'D /'K){q/uj)'^ . The renormalized plasmon frequency 
then follows by solving e{q,Lo) = 1 — VqKeI{{q,Lo) ~ 0. 
For BLG we find that uj^ = ie^erVq- Fig. \Tj^ shows 
renormalized Drude weights V from the full numerical 
calculations described below. 

Optical Conductivity — The full cr^Lo) at arbitrary interac- 
tion strength is obtained by solving the coupled integral 
equations Eqs. ©-([T]) numerically, letting lj Lj + iT~^. 
We can approximately account for screening corrections 
to our first order interaction self-energy by replacing 
the bare Goulomb interaction by its Thomas-Fermi (TF) 
statically screened counterpart. When screening is ne- 
glected, interactions significantly alter ct^lo) in three re- 
spects, as illustrated in Fig. O First, the interband ab- 
sorption threshold is changed dramatically from uj ~ 2eF 
to a; ~ 5£f- This effect is analogous to band gap renor- 
malization in regular semiconductors, with the Fermi 
level playing the role of a gap because of Pauli blocking. 
The threshold shift is equal to 5]^^+ — T,kp- (Fig. [Tj3). 
Second, the value of a{cij) is no longer universal above the 
absorption threshold; instead it shows a decreasing trend 
with to, first reaching above and then dropping below 
the non-interacting value e'^/2h. Third, an absorption 
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peak appears below the threshold. This Mahan exciton 
[2l| feature is a well-understood artifact of our simple 
self-energy approximation. Because of electron scatter- 
ing processes (including Fermi surface fluctuations due 
to intraband electron-hole excitations 2^1 and impurity 
scattering) in the conduction band not captured by our 
self-energy approximation, the Mahan exciton is invari- 
ably unstable. In Fig. [2] we illustrate broadening of the 
Mahan exciton due to disorder. 

The results obtained when we screen the interactions 
in our self-energy expression using a TF approximation 
are shown in grey in Fig. [2l The TF screening wavevec- 
tor for BLG is given by gxF — 4me^/K, a constant inde- 
pendent of electron density. Surprisingly, we find that 
with screening (1) the interband absorption threshold 
shift nearly disappears, (2) RefT(w) ~ e'^/2h above the 
threshold, and (3) the Mahan exciton bound state van- 
ishes. In short, the optical conductivity Recr(aj) behaves 
essentially like that of a non-interacting system. 

To shed light on this result, we observe that the 
TF wavevector qxF — 2.62 x IO^ui^^/k for BLG is 
extremely large compared to all momentum scales of 
electronic transitions and is, in fact, greater than the 
momentum cut-off kc = ^2 77171 for both suspended 
(k = 1) and substrate-mounted (k = 4) bilayers. In 
a TF screening approximation: Vq ~ 27re^/gTF = 
Vo in the regime of interest ep,'-^ < 71, with the 
consequence that in Eqs. ([Sjl- lfTT]) the band renormal- 
ization I]fc+ — Efe_ c>c Vq J d(l)k'kC0s2(f)k'k, and the 
non-equilibrium self-energies i5S^ cx Vq / d(/)fe/feC0S 0fc'fe, 
ST,'^'^ cx Vq J d(j)k'kCOS (j)k'k cos 2 (j)k'k, and iSS"^'^ oc 
Vq J d(/)fe'fcSin0fc/fc sin2(/)fe/fc all vanish. Strong screening 
in BLG restores the optical conductivity essentially to 
its non-interacting value. This remarkable result is pecu- 
liar to BLG, since its double-chirality gives rise to spinors 
with s-wave and d-wave components (rather than s-wave 
and p-wave as in SLG), which do not couple to p-wave 
optical dipole transitions through an s-wave short-range 
interaction. 

Finally we comment on the experimental implications 
of our fin ding s. We interpret the weak experimental 
, [nl, m, [1^ absorption threshold features as evidence 
for short-range screened e-e interactions. We recognize 
that the static screening we use could overstate the re- 
duction in interaction range and that interaction effects 
are likely to persist to some degree, especially in sus- 
pended bilayers for which the dielectric-environment por- 
tion of the screening is absent. Interaction effect could 
be identified experimentally via the Wp renormalizations 
we predict, for example using electron energy loss spec- 
troscopy studies of suspended samples. Interaction ef- 
fects might also be more pronounced in transitions be- 
tween the bands near the Fermi energy - which are in- 
cluded in the massive chiral fermion model we employ 
- and the remote bands located approximately 71 away 
from the Fermi energy. Our theoretical results follow 



principally from pseudospin chirality, and will not be 
strongly influenced by an external potential which opens 
up a gap at the Fermi level of an undoped system. We 
therefore can expect on the basis of current results, that 
excitonic binding energies in these gapped systems will 
be suppressed. 

In conclusion, we have developed a theory for the e- 
e interaction effects on the optical conductivity (t(uj) of 
doped bilayer graphene. We find a novel coupling effect 
which couples the Drude and interband response of the 
optical conductivity, and an accompanying renormaliza- 
tion of the leading-order plasmon frequency. We also find 
that screening dramatically suppresses band renormaliza- 
tion and excitonic effects, restoring a{oj) very close to the 
universal value /2h above the absorption threshold. 
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